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ABSTRACT 

A model is proposed for the formation of water ice mantles on grains in interstellar 
clouds. This occm'S by direct accretion of monomers from the gas, be they formed 
by gas or surface reactions. The formation of the first monolayer requires a minimum 
extinction of interstellar radiation, sufficient to lower the grain temperature to the 
point where thermal evaporation of monomers is just offset by monomer accretion 
from the gas. This threshold is mainly determined by the adsorption energy of water 
molecules on the grain material; for hydrocarbon material, chemical simulation places 
this energy between 0.5 and 2 kcal/mole, which sets the ("true") visible exctinction 
threshold at a few magnitudes. However, realistic distributions of matter in a cloud 
will usually add to this an unrelated amount of cloud core extinction, which can 
explain the large dispersion of observed ("apparent") thresholds. Once the threshold 
is crossed, all available water molecules in the gas are quickly adsorbed, because the 
grain cools down and the adsorption energy on ice is higher than on bare grain. The 
relative thickness of the mantle, and, hence, the slope of Ts^Ay) depend only on the 
available water vapour, which is a small fraction of the oxygen abundance. Chemical 
simulation was also used to determine the adsorption sites and energies of O and OH 
on hydrocarbons, and study the dynamics of formation of water molecules by surface 
reactions with gaseous H atoms, as well as their chances to stick in situ. 

Key words: ISM: dust, ISM:kinematics and dynamics, molecular processes 



1 INTRODUCTION 

The present work addresses three questions: 

1) Why does ice form in clouds only when protected 
from interstellar light by a significant amount of extinction 
7 

2) Why is the optical depth of the ice feature at 3 fi, 
Unearly related to extinction beyond formation threshold ? 

3) Is it possible to quantitatively predict this behaviour 
on the basis of known physics ? 

These questions were addressed more than once in the 
past (see Whittet et al. II2001II . Whittet 1200211 1. However 
the answers given did not always concur. They went roughly 
along two different lines. Initially, van de Hulst 1194911 pro- 
posed that ice formation occurred upon chemical combina- 
tion of O and H atoms on grain surfaces. Later, Jones and 
Williams 1198411 . for instance, developed this model based on 
educated guesses of sticking coefficients of atoms on grains. 

In the meantime, it was also suggested that ice formed 
directly by water vapour deposition. But this model was 
countered by theoretical models of gas chemistry in clouds 
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(Herbst and Leung 119861 1. which pointed to low concen- 
trations (~ 2 10~*) of water vapour in clouds of interest. 
However, addressing the threshold problem, Williams et al. 
1199211 later did not exclude the possibility that this direct 
deposition indeed contributed. Besides, it should be clear 
that this scenario does not in any way imply that all water 
molecules are formed in the primordial gas. 

The present work started in an effort to confirm the first 
condensation scenario by ascertaining quantitatively the es- 
timates which were previously made of the various surface 
parameters involved. Along the way, it was found that wa- 
ter formed on a surface could not stick to it, because of its 
very large initial energy content (~5 eV). This made it dif- 
ficult to pursue the first scenario. Besides, it was also found 
that depositing the energy of one 3-^ photon (~0.4 eV) in 
an adsorbed water molecule, as proposed by Williams et 
al. 119921 1. is not enough to exctract it from the grain sur- 
face, making it necessary to seek another explanation for the 
threshold extinction. Hence the elaboration of a new model, 
based on the second scenario. 

In the following, we reverse this chronological order, and 
start with a detailed description of the proposed model (Sec. 
2), a cursory description of which follows. The first layer of 
ice-to-be develops when the abundance of water molecules in 
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the ambient gas exceeds the vapour pressure corresponding 
to the grain temperature (supercooling). The latter, T, is de- 
termined by equating radiative powers absorbed and emitted 
by the grain: it is, therefore, dependent on the visual magni- 
tude, Av, characteristic of the extinction which shields cloud 
grains from outside radiation. As to the vapour pressure of 
water at this temperature, it depends essentially on T and 
on the heat of adsorption, Bo, of a water molecule on the 
grain surface considered. at threshold is thus indirectly 
linked to Eo- 

Once the first layer is completed, further mantle growth 
is limited only by the total water vapour available in the am- 
bient gas (Sec. 3). This is because, as ice is deposited, the 
grain temperature falls while Eq rises, so that sublimation is 
exponentially suppressed. For a given grain radius, the opti- 
cal depth , ra, at the peak of the ice feature, is then shown to 
be proportional to the increase of extinction over its thresh- 
old value. A relation is given between the corresponding co- 
efficient of proportionality and the relative mantle thickness. 

Assessing quantitatively this scenario (Sec. 4) requires 
the knowledge of the optical and surface properties of grain 
core and ice mantle, as well as heats of adsorption of wa- 
ter molecules on core and mantle. Laboratory data relative 
to the latter is inadequate: most studies are dedicated to 
adsorption of atoms on pure and regular surfaces of tech- 
nical interest, at temperatures and pressures much higher 
than relevant here. However, the availability of new, pow- 
erful, software for chemical modelling has opened other av- 
enues from which to seek the required answers. Adsorption 
(see Steele II2002L Al-Halabi et al. Il2004t 'l. surface dynamics 
(see Kolesnikov et al. (I2UU4II . Hasnaoui et al. 121)1)511 crystal 
growth (see Henkelman et al. Il2003ll ')and chemical reactions 
between gas and solid surfaces (see Sadlej et al. 1199511 ) are 
now abundantly treated with such software. 

The present work also resorts to numerical modeling 
to find adsorption sites and energies for H, O and H2O 
on aromatic and aliphatic hydrocarbon surfaces typical of 
C-rich grains. From this data, and using the condensation 
model proposed above, a threshold extinction can be de- 
duced which is in rough agreemant with observations. 



2 DEPOSITION OF THE FIRST H2O LAYER 

Consider a molecule adsorbed on a uniform flat surface of 
a different material, at temperature T. The probability per 
unit time for it to desorb is 
p = vexp{-Eo/kBT) (s-i) 

where u is the attempt frequency, usually taken to be 
the vibrational frequency of the molecule in its potential 
well, Eq its binding energy or potential depth (also en- 
thalpy of sorption on, or activation energy of desorption 
from, the core surface) and fcs, Boltzmann constant (see 
Atkins II1998II . Somorjai 1199411 1. Assume, ideally, that there 
are s identical adsorption sites per unit area of the surface, 
a fraction 9 of which is occupied by one adsorbate molecule 
each. The desorption rate is then given by 

R = sevexp{-Eo/kBT) (molcm'^s"^) 

If the molecule number density in the atmosphere above 
is fnn, where / is the molecule relative abundance and 
nn, the atomic hydrogen density, and if we neglect for 
the moment the probability of a newly arriving molecule 



being adsorbed over a previously adsorbed molecule (see 
Discussion below), then the deposition or accretion rate is 
fnHV{l — 9)/3, where V is the velocity of the gas molecules, 
and the condition for dynamical equilibrium at coverage 6 
of the surface is 

l/3fnHV{l ~e) ^ s9uexp{-Eo/kBT). (1) 

The mantle formation can be considered to have started 
when ~ 0.5. However, for the first monolayer to be com- 
pleted in a reasonably short period, the desorption rate 
should be, at most, say, one tenth of the deposition rate. 
By definition, this sets the supersaturation ratio at 10. The 
threshold temperature is therefore obtained from eq.(l) by 
setting 6 = 0.5 after multiplying the l.h.s. by 1/10: 

Tthr = ^ , (2) 

where Rb = 1-98 cal/mole °K and Kav = ln{?,Osv / fnuV), 
where the subscript "av" stands for "available" density of 
gas molecules. For the relevant values of Eq, and because 
of the exponential, T is only very weakly sensitive to rela- 
tively strong variations in the components of Kav ■ In zeroth 
approximation, therefore, the poorely known actual abun- 
dance, /, can be replaced by the total oxygen abundance in 
the gas, although the available amount for ice condensation 
is known to be much smaller, because of oxygen depletion 
in CO, OH and grains. Eq, s and v depend on the nature of 
the grain surface and will be estimated in Sec. 4. 

Now, the temperature T of a grain inside a molecular 
cloud is essentially determined by the amount of radiative 
energy falling upon the grain, and this in turn depends on 
the extinction which screens the grain from the outside radi- 
ation. Both relations pave the way to the threshold extinc- 
tion for the onset of ice formation. 

Let the radiative flux at the grain location be F = GFo, 
where _Fo is the interstellar radiation fleld (10~^ Wcm~^) 
and G, the extiction coefficient, which is approximated by 
G = ea;p(-^„/1.08) (see Mezger et al. JT^Hll). 

While the grain, supposed to be spherical, of radius a, 
is still bare, its equilibrium temperature is determined by 
equating absorbed and emitted radiative powers or, approx- 
imately 

GFoTTa?oiy2a = ATTa^2airaasBT'^ , (3) 

where asB is Stefan-Boltzmann constant (5.7 10^^ 
ergcm"^ ° K~*); a„ and air are, respectively, the absorption 
coefficient of the grain material in the visible range and in 
the infrared range corresponding to grain emission at tem- 
perature T. This approximation assumes that the optical 
depth of the grain is that of a cube of side 2a, noting that 
both sides of the equation differ from this by the same factor. 
Now, air can be approximated by its value at the wavelength 
of peak black-body emission at T, i.e. Xmax ~ 3000/T (in 
/i and deg K). Moreover, it is generally found that, in the 
far infrared, which is the emission range of interest here, 
air oc 1/A; that is also the case in the interstellar medium 
(see Dupac et al. ()2()()3^ '). Thus, within the present degree 
of approximation, 

A™„^ - 3000 ' 

where the subscript refers to an arbitrary wavelength 
within the infrared emission range of interest. As a result, 
eq.(3) reduces to 
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Figure 1. Adopted relation between the temperature T of a grain 
and the visual extinction, Av, in front of the grain location in the 
cloud; for 3 sets of values of the grain optical constants, in eq. 4, 
corresponding to 0=19/10^^/'' (lower curve), 19 (intermediate) 
and 19 IQl/^ (upper). 



T = CG^^\ with c = io(i^y/\ 



(4) 



Using the definition of G above, 
T = Cexp{-^J. 

Reasonable values for the grain optical properties are, 
for instance, 

a„ = 2 10* cm"\ Ao = lO^i, ao = 100 cm"\ 

With these values inserted, eq.(4) gives C = 19, so that 



(5) 



a relation plotted in Fig. 1. The latter also shows the impact 
of uncertainty in the optical constants. 

Finally, substituting in eq.(2), the threshold extinction 



Figure 2. The threshold visual extinction as a function of 
the adsorption energy, according to eq.(6), for 5 sets of values 
of adsorption parameters, corresponding respectively to Ka^ = 
Kav,mtidian + NlnW, with N= -2, -1, 0, 1 and 2, from bottom to 
top; N=0 corresponds to the median values defined after eq.(6). 
In all cases, the optical parameters are the same: C = 19. 
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The physics of the condensation process is contained in Kav 
and Eo . These depend on the particulars of the interactions 
of water molecules with the given dust, and are not directly 
available from the literature. Their numerical estimation will 
be taken up in Sec. 4. Here, in order to get a feeling of the 
relevance of this model, we shall study the dependence of 
Av,thr on Eo, considered as a free variable, and on Kav, 
considered as a parameter determined for a few sets of values 
of the physical quantities entering in eq.(6). 

Since there are 5 such quantities, we chose a reasonable 
median average value for each and then add to Kav a number 
of times ±lnlO, thus allowing for large variations of any of 
the 5 quantities. For the median values, we take 

s = 10^^^/ = 510",n// = 10^ / = 7 10'^, 

to get Kav,mcdian = 58.3 and, hence,with C = 19, Fig. 
2. Clearly, the particular choice of these values has no great 
impact on the behaviour of Av,thr as a function of the ad- 
sorption energy. More weight is carried by the optical pa- 
rameters entering in the constant C (eq. 4), although this 
is mitigated by the exponent 1/5: Fig. 3 shows the effect of 
muhiplying the median C by 10"^^^ or 10^^^. 

The horizontal line Av.thr ~ 3, in Fig. 3, corresponds to 
the Taurus cloud as observed by Whittet et al. 11988 ) . Obvi- 



(cal/mole) 

Figure 3. The threshold visual extinction as a function of the 
adsorption energy, according to eq. (6), with the same, median, ad- 
sorption parameters for the 3 curves, but different optical parme- 
ters: C = 19 * 10"^/^, 19 and 19 * 10^5^ from bottom to top. 



ously, such a value can be predicted by the present model, on 
the basis of reasonable or available values of the parameters 
involved. 

However, a wide range of threshold values , from 3 to 30, 
has also been observed (see Williams et al. II19921 '1. which 
prompted the latter authors to propose a model based on 
water desorption upon absorption of resonant infrared pho- 
tons (3 jj.) from the iinterstellar medium. This is discussed in 
Sec. 4, but, in the framework of the present model, a different 
explanation is required. This can be sought in the definition 
of the threshold. In the above, it was tacitly assumed that 
the density distribution within the cloud was single peaked 
and asymmetric (like curve a. Fig. 4) , with Av corresponding 
to the matter between and L/2, distances being counted 
from the star, along the sightline, and ice first appearing at 
the peak, at L/2. It seems much more realistic to envision a 
nearly symmetric distribution (Fig. 4b), in which case, the 



4 R. Papoular 



T.IA, 





11 V \! 

c 


\i \i ] 






a : 






< 1 






\ > 


star / 






\ observer 



L/2 



e/a 



Figure 4. Possible sciicmatic distributions of optically absorbing 
matter along a sight line through a molecular cloud. 



Figure 5. The relative mantle thickness e/a as a function of 
T3/AAv). 



"apparent" threshold is 2Ay^thr- Even this cannot always be 
the case: more probably (see Ward-Thompson et al. ll2UU3ll l. 
the matter in the inner cloud will contribute, in the eye of 
the observer, additional extinction which bears no relation 
to 2Av^thr, which is contibuted by the two limbs (Fig. 4c). 
We therefore expect Av,thr to set only a lower limit to the 
measured, "apparent" threshold. 

Although the other, incompletely known, parameters 
certainly contribute some dispersion of observed thresholds, 
they can hardly explain a wide dispersion among clouds 
which are supposed to harbour a ubiquitous mixture of 
known dust types. 



3 STEADY-STATE ICE MANTLE 

Whittet et al. II1988II showed that, along various lines of sight 
through the Taurus cloud, the optical depth at the peak of 
the ice feature, ra, was a linear function of the corresponding 
extinction: 

T3 = q{A.u - A„,th,)> q=0.072±0.002. 

Such a regular relation seems to imply that a) the grains 
have accreted a constant relative amount of ice, b) this 
amount is the same throughout the cloud. This strongly 
suggests that all water molecules formed through gaseous 
or surface reactions in the cloud are quickly deposited on 
the grains in the inner cloud, which are shielded by A^^thr 
or a still higher extinction, and are therefore cooler still. For 
a given cloud, the number of available water molecules is 
essentially the complement of the number of oxygen atoms 
which are included in CO, OH and grains; it is thus a frac- 
tion of the O abundance. Let us estimate this fraction on 
the basis of observations. For grains of radius a, uniformly 
covered by a mantle of thickness e, the abundance of avail- 
able water molecules, relative to that of H atoms, assumed 
all to be deposited in the mantle, is 

_ 47v{a + e)'^edicengr _ 1 Md dice e , 6 2 /„n 

/ r' ^ J ^ V "!"/' V'/ 

nH^onn o Mh dgr a a 

where Md/Mn is the dust-to-gas mass ratio, dice and dgr 
the mass densities of ice and grain materials respectively. 



On the other hand, for the same geometry and taking 
Nh = 2 lO'^^AAv, where AAv = A^ - A^.thr, 



T3 _ ■R{a + eY .2eaice-ngr-L 



5 10 



-3 Md aice 
Mh dgr 



-(l+-)'-(8) 
a a 



Combining the last two relations and taking Md/Mn ~ 
1/150, dice = 1 and aice = 2.6 10'' cm~' (Leger et al. 
II1983II '). one finally obtains 



/ = 1.3 10 



-3 T3 

Ay 



(9) 



For q=0.072 (Taurus cloud, see above), / = 9.4 lO"'', com- 
fortably smaller than the total O abundance. More gener- 
ally, this relation can be used for a better approximation 
than that which was adopted to plot Fig. 2. 

Eq. (8) also gives e/a as a function of q, a relation plot- 
ted in Fig. 5, with aice = 2.6 lO'^cm"^ and dgr — 2. For 
q = 0.072, e/a = 0.65. 

The assumption that all available water molecules are 
promptly deposited as a mantle on grains in the cloud core 
requires some justification in view of the possible obstacle 
of surface tension. The tension over the curved surface of 
the mantle generates inside an additional pressure which 
enhances sublimation from the surface. This effect is greatest 
at the onset of mantle formation, when the outer radius is 
smallest, and may thus stop deposition altogether if the core 
radius is too small. The critical radius for this to occur is 
deduced from Lord Kelvin's expression for the ratio, Sk, of 
water vapour pressure over a sherical grain of ice to that over 
a plane of the same material. If partial vapour pressures are 
expressed in terms of hydrogen density. 



^ -' /(t, 00, ice) 



exp 



2-f{ice)mH^Q 
d{ice)kBTa ' 



where 'y{ice) is the surface tension, or capillary con- 
stant, of ice. By definition, this is the energy required to 
create a unit surface area of free ice surface, by breaking 
all bonds along a plane drawn through the bulk material. 
This can be approximated by 7(ice) = sEo{ice)/2, where 
Eo{ice) is the adsorption energy of the vapour monomer on 
the ice surface (which is close to the bond energy between 
monomers in bulk ice) , and s, the surface density of adsorp- 
tion sites. Whence, 
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smH20 i^o{^ce;cal niole) 
S,i^ce) = expi^^ ^ ). (10) 

Now, the first fraction on tlie r.li.s. is notliing but tiie 
ratio of surface to volume densities of water monomers in ice, 
i.e., the average distance, 5, between monomers, assumed to 
be the same on the surface and in bulk. This gives finally 

where Eq is in cal/mole. The mantle will grow only if 
the actual pressure of water vapour in the atmosphere, /, 
exceeds f{T, a, ice). Now, in Sec. 2, we defined threshold by 
the condition 

f = Wf{T,Eo{ads)), 

where "ads" serves only as a reminder that the bond 
energy to be used here is that of a monomer on a bare grain 
surface as opposed to Eo{ice), which characterizes bulk ice. 
The condition for the ice mantle to grow can now be written 

Wf{T,Eo{ads)) > Sk{ice)f{T,oo,ice), 

where the f's are given by eq.(l) and depend on the re- 
spective -Eo's but not on a. Taking logarithms of both sides, 

2 3 _|_ (-Eo(iee)-Eo(ads) y S / Eo(ice) 
2T — 2a ^ T 

The adsorption bond energy will be shown below (Sec. 
4) to range between 0.5 and 2 kcal/mole, while Eo{ice) is at 
least 5 times larger (because of the strong polarity of water), 
so that the term 2.3 can be neglected in comparison. Then 
the energies in the other terms approximately factor out, so 
that the condition for mantle growth reduces roughly to 

a>S. 

This growth condition is trivial, and made so because 
the adsorption energy of the monomer on the bare grain is 
so much smaller than the bond energy in the bulk, so that 
the condition for first layer deposition automatically satisfies 
the surface tension constraint. 

It is also necessary to check that the time required for 
mantle deposition is only a fraction of the cloud age. From 
eq.(l). Sec. 2, the time for deposition of a monolayer is 

In the Taurus case, / was found above to be ~ 10~*. 
Assuming uh ~ cm"'^ and T = 10 K at threshold, tiL 
is found to be ~ 10^ y. 

As for the complete mantle deposition time, note that 
the water vapour flux at threshold is /nnV/S; on the other 
hand, in general, the bare grain optical cross-section is 
~ 10~^^ cm^, so the grain surface available for deposition 
is a ~ 4 10~ ^cm per free hydrogen atom (Spitzer (|1978jl '). 
Thus, the total amount of water ultimately deposited as ice 
is fnnVa/S per second per H atom, while the total number 
of water molecules assumed to be available in the atmo- 
sphere is / per H atom. The time required for deposition is, 
therefore, 

tmantle ~ ^nVa (^)' 

For the same uh and t as assumed above, this yields 
~ 3 10^ y. Note that this time is independent of grain ra- 
dius, mantle thickness or water vapour abundance in the 
atmosphere. 

Thus, both deposition times are reasonably short 
enough for the model to apply. This model can be used for 
predictive purposes. In the particular case of the Taurus 
cloud, for instance, the observations of Whittet at al. II1988II 
suggest that the "true" threshold is, at most, 3.3 mag. Fig- 
ure 3, here, indicates that this range corresponds roughly 
to 1000 < Eo < 2000 cal/mole, which may help charac- 



terize the dust on which ice condensed. Also, since no ice is 
positively observed in the diffuse interstellar medium, where 
Av = 0, the same figure predicts that no dust material has 
Eo > 3 kcal/mole. Vice versa, if that dust is fully character- 
ized beforehand. Fig. 3 will predict the "true" threshold for 
a given Eq. Either way, prediction requires prior knowledge 
of the adsorption bond energy of water on the candidate 
dust material. Unfortunately, very little seems to have been 
published so far on the materials of interest to astronomy, 
which brings us to the subject of Sec. 4. 



4 ATOMIC, MOLECULAR AND SURFACE 
PROPERTIES 

The discussion above emphasized the need for quantitative 
data on which to base and justify a model. Relevant data 
in the literature are meagre at present, most of the recent 
effort being directed to interactions of atoms with pure, reg- 
ular surfaces of materials of technical interest (see Somorjai 
II1994II . Bortolani et al. (119941 . Masel J1996II '): water and ice 
are certainly not a priority even though an impressive spe- 
cialized review of this case was given some time ago by Thiel 
and Madey II1987II . Masel points to the odd behaviour of 
H2O and CO, which do not follow general trends; he gives 
typical physisorption energies as 2 to 10 kcal/mole. Dormant 
and Adamson (1968) studied the adsorption of nitrogen on 
molecular solids. For ice, in particular, they give Eq between 
1500 and 2500 cal/mole. 

More to our point, Zettlemoyer et al. Ill975t studied 
the adsorption of water on several solid surfaces, including 
carbon blacks, for which they find adsorption energies ~ 10 
kcal/mole. But Suzuki and Churchill 1199011 cited a range 
of 1 to 10 kcal/mole for common gases on sieving carbon. 
Ibupoto and Woods 1200111 measured about 600 cal/mole for 
water on cellulose. 

These widely dispersed results were also generally ob- 
tained at about 100 K and relatively high pressures. They 
can hardly be extrapolated to the special cases of adsorption 
of O, OH and H2O on amorphous ice and hydrocarbons at 
much low temperature and pressure (see further discussion 
in Sec. 5). 

Of course, new experiments will ultimately give the an- 
swers, but numerical chemical simulation is now able to de- 
liver more quickly an estimate of the required quantities 
relevant to the special case at hand. It is used here, instead. 

4.1 The chemical code 

A description of the software package used for simulation 
can be found in Papoular 1200H . The particular code used 
here is AMI (Austin Model 1), as proposed, and optimized 
for carbon-rich molecules, by Dewar et al. (1985 1. This semi- 
empirical method combines a rigorous quantum-mechanical 
formalism with empirical parameters obtained from com- 
parison with experimental results. It computes approximate 
solutions of Schroedinger's equation, and uses experimen- 
tal data only when the Q.M. calculations are too difficult 
or too lengthy. This makes it more accurate than poor ab 
initio methods, and faster than any of the latter. 

The code is used in the first place to create a plau- 
sible model grain: once the size, composition and general 
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structure have been chosen, the "optimization" procedure 
determines the correct relative positions of the atoms by 
minimizing the total potential energy. In order to study the 
interaction of a given species (atom or molecule) with the 
grain surface or with species previously adsorbed on it, the 
new species is created at some distance from the grain and 
given some velocity towards it. The molecular dynamics is 
then launched. 

The molecular dynamics relies on the Born- 
Oppenheimer approximation to determine the motions 
of atoms under nuclear and electronic forces due to their 
environment. At every step, all the system parameters 
are memorized as snapshots so that, after completion of 
the run, a movie of the reactions can be viewed on the 
screen, and the trajectory of any atom followed all along. 
This makes for a better understanding of the details of 
mechanisms and outcomes. Note that the dynamics of 
bond dissociations and formation can only be simulated by 
using Unrestricted Hartree-Fock (UHF) wave functions in 
the Q.M. part of the calculation (see Szabo and Ostlund 
Ill989t '). The elementary calculation step was set at 10~^^ 
s. 

Semi-empirical simulation methods do not account for 
purely quantum-mechanical phenomena like tunneling or 
zero-point energy. In the present context of relatively high 
grain temperature, the former is negligible as regards surface 
diffusion, and the latter hardly affects the probability and 
energetics of the physico-chemical reactions involved. The 
procedure just described was used by Papoular 12UU5II to 
study the interactions of atomic H with hydrocarbon grains. 
In the course of that work, the adsorption potential wells of 
the model grains were determined in position and depth. 
The present work made use of these results as a guide for 
the adsorption of O, OH and H2O. 

4.2 Adsorption sites and energies of H2O 

This study is restricted to carbonaceous IS (InterStellar) 
dust. IR spectra of such dust reveals the presence of both 
aromatic and aliphatic C-C and C-H bonds. In order to save 
computation time, we have modeled several small structures 
presenting the most typical functional groups. Adsorption 
sites can be found by tossing the water molecule towards 
one of these targets, with some initial velocity. If this veloc- 
ity is small enough, the projectile will dwell on the target 
periphery and wander over it, progressively loosing energy 
to the target, until it settles at the bottom of one of the 
possible adsorption wells created by the spatial variations of 
the van de Waals potential of the structure. 

A quicker procedure is based on the use of the previous 
results on H adsorption: since the sites roughly coincide, the 
water molecule is initially positioned in one of the known 
sites of H adsorption, all velocities being set to zero, and 
the potential optimization of absorbant and adsorbate to- 
gether is then launched; note the (possibly) new position of 
the adsorbate, and let Emin be the minimum potential thus 
found; now, extract the water molecule and read the new 
potential of the structure, say E grain- Then, the adsorption 
energy can be deduced to be 

Eq = Egrain + EH2O ~ Emin, 

where Eh^o is the internal potential energy of a water 
molecule, and is computed by the code to be 223 kcal/mole. 




Figure 6. Main adsorption sites and energies near typical liydro- 
carbon surface functional groups. Aromatic CH groups are recog- 
nized, as well as aliphatic CH, CH2 and CH3 groups, in the con- 
ventional representation. In each case, the independent V-shaped 
molecule is water, trapped in the adsorption site, and its oxygen 
atom is linked to the nearest H atom or atoms on the "grain" by 
dotted lines whose lengths are indicated by near-by numbers in 
angstroms. Below each system appears its ranking number, fol- 
lowed by the corresponding adsorption energy, Eq, in kcal/mole. 
One typical site not represented in this figure is above the center 
of aromatic rings, at ~ 3.5 A from the 6 carbon atoms, with the 
lowest adsorption energy: 0.5 kcal/mole. (1 kcal/mole=43.5 meV) 

Egrain usually docs not differ much from the potential en- 
ergy of the target alone, because van der Waals forces are 
so weak relative to the chemical binding energies. 

The initial temperature of each system can be controled 
through the chemical code. Here, in view of the very low 
temperatures known to be required for ice formation, the 
initial temperature was uniformly set at K. When the 
water molecule is adsorbed, energy Eq is shared between 
adsorbate and absorbant; because the latter is normally so 
much larger than the water molecule, the consequent rise in 
temperature is negligible and the adsorbat remains locked 
to the surface. 

Figure 6 summarizes the result and shows that the ad- 
sorption energies determined as described fall between 0.5 
and 2 kcal/mole; according to Fig. 3, the proposed model is 
therefore coherent with observations. It also shows that, be- 
cause these energies are so weak, ice can hardly form in the 
diffuse interstellar medium, and needs at least a few visible 
magnitudes of shielding. This result is discussed further in 
Sec. 5. 

It is noteworthy that an average bonding energy be- 
tween monomers in ice can be deduced from the latent 
heat of sublimation of ice in vacuum, 600 cal/g, to be 
18*600*10"^ = 10.8 kcal/mole, much larger than adsorp- 
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tion energy of a monomer on a bare grain, and nearer to 
the results of temperature programmed desorption (see Kay 
et al. I|1989|l . Eraser et al. ("2001)). In order to understand 
this behaviour, successive layers of monomers were added 
on each other to form an ellipsoidal chunk of 27 molecules. 
By minimizing its potential energy at zero temperature, this 
chunk was "quenched" into an amorphous state in which all 
monomers fit tightly, but randomly, together. Starting from 
the outer ones, the monomers were then extracted in turn, 
noting the new potential energy at each step. From this it 
was deduced that the binding energies ranged from 18 to 2 
kcal/mole, or roughly 10 on average. Ice monomers hold to- 
gether more tightly than to the substrate, because hydrogen 
bonds are stronger than van der Waals forces. 

Further quantitative information can be gained from 
the chemical model used in this section. Thus, based on the 
average distance of adsorption sites from the outer atoms of 
the grain ( Papoular (2005)), we estimate the average dis- 
tance between these sites to be about 3 A . This justifies 
setting s = 10^^ cm~^ for the surface density of states in eq. 
(1). Moreover, once a water molecule has settled at the bot- 
tom of one adsorption potential well, it can be displaced by, 
say, ~ lA and then left to oscillate about the bottom while 
monitoring the dynamics all along by means of the chem- 
ical software. This gives us the attempt frequency, v, used 
also in eq. (1), about 5 10^^ s~^, on average. This value falls 
in the range typical of (the relatively weak) physical bonds 
,~ 10^^ s~^, as opposed to ~ 10^'' s~^, which applies to (the 
stiffer) chemical or even H-bonds ( see Atkins lH^^ 'l. 

4.3 Ice formation on grains 

Chemical simulation also helps discussing the possibility of 
formation of ice on grains by accretion and conversion of O 
atoms into water, followed by combination with two succes- 
sive gaseous H atoms impinging on the solid surface. For this 
purpose, the procedure used previously for H adsorption, 
and above for water adsorption, was applied to O atoms. 
Here, for a more realistic model, a larger structure (18 C 
and 28 H atoms), similar to that used by Papoular((2flQy ), 
and including the functional groups of Fig. 6, was used in- 
stead of the individual functional groups. Adsorption sites 
were found along CH functional groups at ~ 2.2^4 from 
the H atom, with energy 0.6 kcal/mole. This is enough 
to keep the adsorbate on the surface at the low tempera- 
tures of interest, i.e. the sticking probability can be taken 
to be 1. Next, a free H atom was directed roughly towards 
an adsorbed O atom. After some wandering it reacted with 
the latter to form an OH radical which was immediately ex- 
pelled, with a high velocity (1.5 10^ cms~^) and strong in- 
ternal oscillations (70 kcal/mole), both tapped from the high 
chemical energy released in the process (111 kcal/mole). Ex- 
pulsion occured within 1 period of internal oscillation of the 
radical (^ lO"'^* s), which is much shorter than the relax- 
ation time of the grain (10~^^ — 10~^^ s). 

The same procedure was applied to the second step of 
this scenario, viz. accretion of a gaseous H atom by an ad- 
sorbed OH radical, to form water. Here, by contrast with 
OH formed in situ, the OH radical is first independently op- 
timized in potential energy, and deposited with zero atomic 
velocities near an adsorption site. The whole system is then 
optimized to find the adsorption energy, ~ 1 kcal/mole. As 



before, a free H atom was then directed towards the ad- 
sorbed OH radical, and ultimately formed with it a new 
water molecule. The latter was expelled from the surface 
with high internal and translational energies, but, this time, 
only after ~ 2 ps and a large number of internal oscillations; 
by then, the grain had relaxed. 

Both the expulsions of newly formed OH and H2O can 
be understood by monitoring the motions of their atoms as 
well as of the nearest H atoms of the target grain. When a 
new OH bond is formed, the H atom oscillates with a large 
amplitude and, in so doing, encounters surface H atoms. 
At each encounter, part of its internal energy is lost to the 
grain and another part is transformed into external kinetic 
energy of the new species. When the latter becomes about 
equal to the adsorption energy of the species on the grain, 
the species escapes. Expulsion is favored by low adsorption 
energies and high formation energy of the new species. The 
grain size does not matter as long as it is much larger than 
the said species. 

Thus, although water molecules may form on grains ac- 
cording to van de Hulst suggestion, neither the intermediate 
OH, nor the ultimate H2O can stay in situ. Direct accre- 
tion of water from the gas is required to start ice mantle 
formation. 

4.4 Photodesorption of water 

Williams et al. II1992II proposed that, if enough interstellar 
radiation at 3 is allowed on adsorbed water molecules, 
they can be desorbed because this wavelength is resonant 
with their natural frequency, and can therefore impart much 
energy to them; and this would explain the existence of a 
threshold for ice formation. In order to test this hypothe- 
sis, a water molecule was deposited with zero velocity in 
an adsorption site; one of its OH bonds was then extended 
from its equilibrium length (0.96 A ), to 1.2 A , which gave 
the molecule an excess potential energy of 10 kcal/mole (the 
energy of a 3/x photon). Strong oscillations set in, but the 
molecule was not desorbed. The initial excess energy was 
then increased by steps till the molecule was observed to 
desorb. This did not happen before the excess reached near 
the OH bond energy (111 kcal/mole), in agreement with 
the finding above that OH or H2O formed on the surface by 
combination of O or OH with a free H atom, is immediately 
desorbed. Thus, a single infrared photon can hardly desorb 
an adsorbed monomer. 



5 DISCUSSION 

Section 2 addresses only the very beginning of mantle forma- 
tion, when water has not yet filled all the available adsorp- 
tion sites on the bare grain surface. Equation 1 assumes that, 
in this sub-monolayer regime, the desorption rate is propor- 
tional to the surface coverage by water molecules. In general, 
this "first-order" model is considered to be legitimate in as 
much as it is applied to the initial stage of adsorption, under 
very low temperature and vapour pressure, when the grain 
is nearly bare and the molecules hit it one at a time and stay 
put wherever they are adsorbed. Under such circumstances, 
before a large fraction of the grain surface has been covered 
with a monolayer of adsorbate, a newly arriving molecule 
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stands only a minute chance of falling precisely on a previ- 
ously adsorbed one, so that adsorption and desorption are 
essentially independent processes, both being, however, con- 
troled by interactions of a single adsorbate molecule with 
the bare surface. This is the basis of Langmuir's adsorption 
isotherms, which is found to apply most generally, especially 
to non-metallic substrates (see Masel ( 1996 1, Atkins ( 1998}). 

Later on, incoming molecules pile up upon one another, 
forming a multilayered structure. For water and alcohols, 
prone to H-bond formation, the measurements of desorption 
rates from such multilayers over plane substrates generally 
point to plane-like layers, from which desorption proceeds 
"layer- by-layer" , according to a so-called "zero-order" rate 
(e.g. Kay et al. Ill9ij9ll . Fraser et al. II2UU1L for water ice); 
this is to mean that the desorption rate is now independent 
of coverage, unlike the sub-monolayer case. Cases have been 
observed where a fractional-order dependence is observed 
(see Nishimura et al. (19981; Wu et al. ilggg} !, for alcohol 
ices) . This multilayer regime is not relevant to the threshold 
problem but applies to ice mantle formation (Sec. 3): a non- 
zero desorption order (which does not seem to have been 
observed for water) could slow down the growth rate, but 
only slightly, for a given multilayer adsorption energy. 

Nontheless, in order to understand published exper- 
imental results, it may be useful to dwell on multilayer 
structures, which are most often encountered in the wa- 
ter literature. In such structures, adsorbate-adsorbate in- 
teractions dominate over adsorbate-substrate bonding. As a 
consequence, both _Eo and the pre-exponential factor in the 
desorption rate may differ from those of the sub-monolayer 
regime. For hydrogen-rich species like water or alcohol ice, 
hydrogen bonds between adsorbates make for values of Eo 
higher than those typical of the common van der Waals 
bonds responsible for physisorption. Thus, while Eq is of 
order 10 kcal/mol for internal bonds of water and alcohol 
ices (see above and Fraser et al. 1200111 and Nishimura et al. 
H1998ll 'l. it is reduced, for instance, to about 2 kcal/mol for 
liquid water, in which H-bonds do not form. 

In general, when the bond energies of the adsorbate to 
itself and to the substrate are distinct but slightly different, 
both can be deduced by analysis of the two corresponding 
peaks of desorption rate as a function of temperature, in 
the multilayer and sub-monolayer regimes, respectively (see 
Redhead HW^ ). For water on a gold substrate, the des- 
orption curves yield only the water-water bonding charac- 
teristics, because they do not display a well-resolved sub- 
monolayer peak. Kay et al. Ill989t interpreted this result as 
an indication that water binds more strongly to itself than 
to gold, to which it does not bond chemically. Zettlemoyer 
et al. 1197511 also found that carbon blacks were hydropho- 
bic. This hydrophobicity strongly favours tight water-water 
bonding over weak water-substrate bonding, even in the 
sub-monolayer regime, the more so the higher the tempera- 
ture and pressure. As a consequence, the adsorbed molecules 
would tend to bind together and segregate in 3-D islands, 
and desorption would, very early on, convert to a zero-order 
course, indistinguishable from the later course. This exam- 
ple shows the difficulty of determining experimentally ph- 
ysisorption bond energies, which requires detailed and ac- 
curate sub-monolayer measurements at very low tempera- 
ture (see Masel JT^Ml V In a rare, but successful endeavour, 
Andersson et al. 1198411 studied the adsorption of water on 



Cu(lOO) and Pd(lOO) around 10 K, at submonolayer cover- 
ages. Using electron energy loss spectroscopy (EELS), they 
showed that water adsorbed as monomers with the oxygen 
end towards the substrate, and identified a hindered transla- 
tional vibration at 28.5 and 41.5 meV, respectively. This cor- 
responds to adsorption energies of 0.65 and 0.95 kcal/mol. 
They also observed the transition to clustered water beyond 
a temperature of 20 K or a coverage of 0.2. 

Our chemical model for adsorption of water on a hy- 
drocarbon substrate yields bond energies between 0.5 and 
2 kcal/mol (depending on the particular site) and, correl- 
atively, reveals no chemical or H-bonds between water and 
the model substrate (i.e. there is no electron sharing between 
water and substrate). From the known performance of the 
AMI class of semiempirical models, the uncertainty in this 
estimate is a factor 2, similar to the intrinsic dispersion due 
to differences between adsorption sites. 

Although it leaves to be desired, our estimate is quali- 
tatively confirmed by the following facts: 

• the finding by Zettlemayer et al. 119751 that carbon 
blacks are hydrophobic, and by Kay et al. 1198911 that water- 
substrate bonds are much weaker than internal water-ice 
bonds; 

• the measurement of adsorption energies of water on Cu 
and Pd by Andersson et al. 1198411 . giving about 1 kcal/mol; 

• the internal liquid water bond energy derived from its 
latent heat (80 cal/g) or surface tension (73 erg/cm^) is 2 
kcal/mol; CO-CO surface binding energy in pure CO ices 
is about the same (Sandford and AUamandola (1988)); Al- 
Halabi et al. 120041 find a range of 1 to 3.5 kc/mol for CO 
on water ice; 

• it is observed in the sky that mantle formation requires 
cooling inside shielding clouds down to about 10 K, while 
desorption of multilayer ice mantles in the laboratory occurs 
consistently above 130 K; if the internal (chemical) bond 
energy of water ice (~ 10 kcal/mol) were to apply to our 
threshold model in Sec. 2 instead of the physisorption bond 
energy determined in Sec. 4 (0.5 to 2 kcal/mol), fig. 2 shows 
that Ay^thr would always be ~ 0, and ice would condense 
even in the diffuse ISM, neither of which is observed. 

A side effect of the high internal binding energy of wa- 
ter, the attendant formation of 3-D "droplets" and the early 
setting of the zero-order regime, is that the threshold con- 
dition should no longer be defined by eq.(l) with 9^1, but 
with a lower value. All other parameters being equal, this 
would decrease Kav in eq. 2, increase the threshold tem- 
perature and hence, decrease the threshold extinction. For 
instance, Nishimura et al. 1199 8 1 for methanol on alumina, 
and Wu et al. 119931 1 for alcohols on NiO, observed the for- 
mation of "droplets" as early as 6' ~ 0.5. For a carbon sub- 
strate, Andersson et al. p984t find the transition at ~ 0.2 
monolayer. Even if 9 were much smaller, that would not have 
dramatic effects on the threshold, thanks to the logarithmic 
expression of Kav (eq.(6) and Fig. 2). As observed before, 
the desorption order in the following multilayer regime is 
not relevant to the threshold problem. 

Among the other causes of uncertainty in our estimate 
of threshold extinction, a minor one is the value of the at- 
tempt frequency of water physisorbed on the model sur- 
face, ly — 510^^^ s"'^ as compared with 5 lO'^^ s~^ for wa- 
ter H-bonded to its ice. Generally, for various materials, it 
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is considered that typical values are of order 10^^ s ^ for 
physisorption and 10^^ s~^ for chemisorption (see Atkins 
JinSHll, Somorjai l|T?M|l . Masel (IHMjl V the strength of H- 
bonds being generally intermediate between those of physi- 
and chemisorption. Temperature desorption spectra yield 
10^^ s~^ for methanol H-bonded on alumina (Nishimura et 
al. (EHHHJ) and lO" s~^ on NiO (Wu et al. llT^ ). Our esti- 
mates fit well in this picture because a weaker bond results 
in a weaker restoring force and, hence, a lower oscillation 
frequency of the adsorbate in its potential well, i.e. a lower 
attempt frequency. Besides, even if the attempt frequency 
were in error by 1 or 2 orders of magnitude, the conclusions 
would not be affected notably, thanks again to the over- 
whelming weight of the exponential in eq.(l) (see fig. 2). 



6 CONCLUSION 

Chemical simulation supports a model of mantle forma- 
tion by direct accretion of monomers from the gas, be 
they formed by gas or surface reactions. It also supports 
the threshold being defined by the minimum visible extinc- 
tion of interstellar radiation which is necessary to lower the 
grain temperature to the point where thermal evaporation 
of monomers is just offset by monomer accretion from the 
gas. 
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